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With wavefunction-based electronic-structure calculations we determine the Cu d-d excitation 
energies in quasi-one-dimensional spin-chain and ladder copper oxides. A complete set of local 
excitations has been calculated for cuprates with corner-sharing (Sr2Cu03 and SrCu02) and edge- 
sharing (LiVCu04, CuGeOs, LiCu2 02 and Li2Cu02) Cu04 plaquettes, with corner-sharing CuFe 
octahedra (KCuFs), for the ladder system CaCu2 03, and for multiferroic cupric oxide CuO. Our 
data compare well with available results of optical absorption measurements on KCuFs and the 
excitation energies found by resonant inelastic x-ray scattering experiments for CuO. The ab initio 
results we report for the other materials should be helpful for the interpretation of future resonant 
inelastic x-ray scattering experiments on those highly anisotropic compounds. 



I. INTRODUCTION 

Many-body effects within the partially filled 3d shells 
of transition-metal oxides are at the heart of various in- 
triguing phenomena. Much of the interesting physics 
comes from the interplay between electron delocalization 
caused by intersite orbital overlap and strong Coulomb 
interactions within the d shell. The investigation of this 
interplay goes back as far as the 1930's with the work 
of Verwey, Mott, and others."^ Additional degrees of free- 
dom arising from orbital degeneracy, electron-lattice cou- 
plings, and/or ligand p to metal d charge transfer effects 
increase in many cases the complexity of the problem. 
Further, it was found that the multiorbital couplings 
within the 3d shell may be strong even when the orbital 
degeneracy is lifted, e.g., by the electrostatic field set 
up by neighboring ions. For the copper oxide supercon- 
ductors, for example, it is believed that the off-diagonal 
coupling between states of x'^—y'^ and z'^ symmetry is one 
of the elements which determine the precise shape of the 
Fermi surface, see, e.g., [2]-[6| and references therein. Ac- 
curate estimates for the size of the d-level splittings are 
therefore important in designing realistic model Hamil- 
tonians for these systems. 

Traditionally the d-d excitations in transition-metal 
oxides have been investigated by optical spectroscopy. 
First measurements on rocksalt compounds such as NiO 
and CoO were reported in the late 50's.^ Although in 
centro-symmetric cases these transitions are not dipole 
allowed, they do aquire weak intensity in the optical spec- 
tra due to the admixture of lattice vibrations which lift 
the inversion symmetry at the transition-metal site. 

It has been shown recently that highly accurate mea- 
surements for both magnetic and local charge excitations 



can be carried out by resonant inelastic x-ray scattering 
(RIXS).^ As demonstrated by Moretti Sala et al^ for the 
case of d-d transitions in two-dimensional (2D) Cu oxide 
systems, the high energy resolution and the analysis of 
the polarization and scattering geometry dependence al- 
low a quite reliable interpretation as concerns the nature 
of the final excited states. 

The calculation of local excitations within the 
transition- metal 3d shell is an interesting problem. It 
is desirable to use alternative approaches to standard 
density functional theory (DFT) calculations because 
in first place electron correlations are strong and sec- 
ondly, in its original formulation, DFT is only a ground- 
state theory. The ab initio wavefunction-based quantum 
chemical approaches therefore constitute the method of 
choice here. Advanced quantum chemical calculations 
have been recently applied to the study of d-d excita- 
tions in La2Cu04, Sr2Cu02Cl2, and CaCu02. The ab 
initio results^^ turned out to be in good agreement with 
the RIXS data reported by Moretti Sala et al.^ On the 
other hand, recent DFT calculations on La2Cu04 pre- 
dict dx2_y2 to dz2 excitation energies that are 0.5 eV 
lower than in experiment.^ ^ 

Here, we extend the study of 2D cuprates as in Ref. |10] 
to the case of highly anisotropic chain and ladder Cu d^ 
compounds. The motivation for the present investigation 
is twofold. One one hand, we want to check the perfor- 
mance of our quantum chemical computational schem^^ 
in the case of one-dimensional (ID) and quasi-lD Cu d^ 
systems. Quantum chemical calculations on relatively 
small clusters have been earlier performed on both 2D 
and ID cuprates, e.g., the 2D layered materials La2Cu04 
and Sr2Cu02Cl2 and the chain-like system Sr2Cu03.'^ 
However, differences as large as 0.5 eV were found be- 
tween the on-site d-d excitation energies reported for 
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La2Cu04 and Sr2Cu02Cl2 in Ref.TTand our more recent 
results discussed in Ref. Those differences seem to be 
related to the less precise description in Ref. [11] of the 
adjacent Sd-metal and O ions, i.e., the nearest neighbor 
(NN) Cu ions and ligands around the CUO4 plaquette 
at which the d-d excited states are explicitly computed. 
Similar differences are here found for the d-level splittings 
of the chain system Sr2Cu03, which shows that a careful 
analysis of this issue is indeed motivated. 

Secondly, our ab initio data should be helpful for the 
correct interpretation of RIXS and optical spectra in 
these compounds. For highly anisotropic structures, the 
degeneracy of both the t2g and Cg levels is lifted and the 
excitation spectra display a very rich structure. Even in 
the 2D compounds, the sequence of the different excited 
states cannot always be predicted beforehand. For exam- 
ple, due to the different ratios between the in-plane and 
apical Cu-0 bond lengths, the hole state corresponds 
to the lowest crystal- field excited state in La2Cu04, with 
a relative energy of 1.4 eV,'^^and to the highest crystal- 
field excitation in HgBa2Cu04, with a relative energy of 
2.1 eV.II^ The situation should be even more complex in 
ID systems. Reliable ab initio results are therefore de- 
sirable for the chain and ladder cuprates. 



II. COMPUTATIONAL DETAILS 

Our computational approach is based on correlated ab 
initio methods traditionally used in quantum chemical 
studies on molecular systems. For each of the materi- 
als investigated here, the starting point is a restricted 
Hartree-Fock (RHF) calculation with periodic boundary 
conditions. All RHF calculations were performed with 
the CRYSTAL program package. We applied Gaussian- 
type atomic basis sets from the crystal library, i.e., ba- 
sis functions of triple-zeta qualitjEH for Cu, O, and F and 
basis sets of either double-zeta or triple-zeta quality for 
metal ions next to the CUO2 chains (e.g., Li, Ca or V). 
In all computations experimental lattice parameters were 
used, as reported in Refs. USHUl 

Post Hartree-Fock correlation calculations can be car- 
ried out on finite embedded clusters, due to the local 
character of the correlation hole of a (i electron. Yet, the 
orbitals used here are those of the infinite system. Each 
embedded cluster C consists of two distinct regions: an 
active region Ca where the actual correlation treatment 
is performed and a buffer region Cb whose role is to pro- 
vide support for the longer-range tails of Wannier orbitals 
(WO's) centered at sites in Ca- The active region Ca in- 
cludes one reference Cu site, the NN ligands, and the NN 
Cu ions. For the systems addressed in this study, a given 
Cu ion may have four, five or six NN ligands. Those NN 
ligands thus form either L4 plaquettes, L5 pyramids or 
distorted Lq octahedra around a particular Cu site. As 
concerns the Cb region, we include in there each ligand 
coordination cage around the NN Cu sites and all NN 
closed-shell metal ions. In CuGeOa, for example, there 



are 8 Ge^+ NN's. In LiCu202, there are 5 Cu^^ 3d^^ 
and 8 Li^+ NN's. All post-RHF computations were per- 
formed with the MOLPRO quantum chemical software.^^ 
The orbital basis associated with a given cluster C is 
a set of projected WO's: localized WO's associated with 
the RHF bands are first derived with the help of the 
orbital localization module of the crystal package and 
subsequently projected onto the set of Gaussian basis 
functions associated with the atomic sites within the clus- 
ter C. Technicalities concerning this procedure are dis- 
cussed in Refs. 10, 23l and[24l For each particular cluster, 
the RHF data is additionally used to generate an one- 
electron embedding potential that models the crystalline 
environment. This effective potential is constructed with 
the MATROP module of the MOLPRO program by using 
a real-space matrix representation of the self-consistent 
Fock operator from the periodic RHF calculation.l^^I^ All 
necessary RHF data is converted into MOLPRO format 
with the help of an interface program. Although the 
WO's at the atomic sites of C are derived for each of the 
compounds discussed here by periodic RHF calculations 
for the Cu 3d^ electron configuration, the embedding po- 
tentials are obtained by replacing the Cu^+ 3d^ ions by 
closed-shell Zn^+ 3d^^ species. This is a good approxi- 
mation for the farther Sd-metal sites, as the comparison 
between our results and RIXS data shows. ^ The exten- 
sion of this embedding scheme toward the construction 
of open-shell embeddings is an ongoing project in our 
group. 

While the occupied WO's in the buffer region Cb 
are kept frozen, all valence orbitals centered at ligand 
and Cu sites in Ca are further reoptimized^^ in mul- 
ticonfiguration complete-active-space self-consistent-field 
(CASSCF) calculations.^^ In the latter, the ground-state 
wavefunction and the crystal-field excited states at the 
central Cu site are computed by state-averaged multiroot 
optimizations.^^ The Cu d-level splittings are finally ob- 
tained from additional multireference single and double 
configuration-interaction^^ (MRCI) calculations as the 
relative energies of the crystal-field excited states. The 
central Cu 3s, 3p, Sd, NN ligand 2p, and NN half- filled 
Cu 3d orbitals are correlated in MRCI. 

The virtual orbital space in the MRCI calculation can- 
not be presently restricted just to the Ca region. It thus 
includes virtual orbitals in both Ca and C^, which leads 
to very large MRCI expansions. To make the compu- 
tations feasible, we restrict our study to ferromagnetic 
(FM) allignment of the Cu d spins. Even for FM clus- 
ters, the MRCI expansion may include in some cases up 
to ~10^^ Slater determinants. 



III. RESULTS AND DISCUSSIONS 

A. Chains of corner-sharing CUO4 plaquettes 

The ID compounds Sr2Cu03 and SrCu02, built of 
chains of corner-sharing CUO4 plaquettes, display a 
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number of quite unusual properties. On one hand, 
these materials have the largest NN antiferromagnetic 
(AF) exchange integrals in the family of Cu ox- 
ides. Yet, the interchain couplings are very weak.l^^M^ 
They thus constitute ideal systems for studying the mag- 
netic response of ID spin-1/2 antiferromagnets. Sec- 
ondly, high-resolution angle-resolved photoemission ex- 
periments on these compounds have for the first time 
revealed the realization of spin-charge separation in ID 
electron systems.-^ 

We list in Table [l] results for the d-level splittings 
of Sr2Cu03 and SrCu02. As concerns the crystallo- 
graphic details, a major difference between these two 
compounds is the presence of a network of single Cu02 
chains in Sr2Cu03,^^ see Fig. 1, whereas SrCu02 displays 
a double-chain structure, with CUO4 plaquettes on adja- 
cent chains sharing common edges. The clusters we use 
in our CASSCF and MRCI calculations include there- 
fore two NN plaquettes in the case of Sr2Cu03 and four 
adjacent plaquettes for SrCu02. The reference system 
is chosen such that the x-axis is along the Cu02 chains 
and z is perpendicular onto the CUO4 plaquettes. In 
this framework, the ground-state hole orbital has —y'^ 
symmetry. The data in Table [l| show that the (i-level 
splittings in Sr2Cu03 and SrCu02 are very similar. In 
both cases, the relative energy of each particular excited 
state is not very different from the corresponding value 
in the Li2Cu02 compound, see Sec. III.C. As in Sr2Cu03 
and SrCu02, there are no apical ligands in Li2Cu02 ei- 
ther. Results for the d-level electronic structure of the 
latter compound are listed in Table [TTI| 

Inclusion of single and double excitations on top of the 
CASSCF wavefunctions brings a nearly uniform upward 
shift of 0.1 to 0.2 eV for the d-d transitions. By adding 
Davidson corrections,'^ the relative energies of the ex- 
cited states further increase by 0.05 to 0.1 eV. The MRCI 
values listed in Table [TTll include such Davidson correc- 
tion terms. This relative stabilization of the electronic 
ground state with respect to the higher lying states in 
MRCI is mainly related to ligand 2p to metal 3(i charge- 




FIG. 1. Chain of corner-sharing Cu04 plaquettes in 
Sr2Cu03. The ground-state hole orbital is (i^2_^2. 



TABLE I. CASSCF/MRCI d-d excitation energies for corner- 
sharing chains of Cu04 plaquettes in Sr2Cu03 and SrCu02 
(eV). The MRCI values include Davidson corrections."^^ 



Hole orbital 


SrsCuOa 


SrCuOs 




CASSCF/MRCI 


CASSCF/MRCI 


2 2 
X -y 








xy 


1.20/1.50 


1.26/1.55 


yz 


1.83/2.14 


1.77/2.03 


xz 


1.90/2.24 


1.88/2.19 


z" 


2.16/2.55 


2.08/2.44 



transfer correlation effects which are more important for 
the in-plane d^-^-y^ hole orbital having a- type overlap 
with the NN O 2p functions. 

As concerns the comparison with other theoretical in- 
vestigations, there are differences of 0.3 eV or more be- 
tween the d-d excitation energies listed for Sr2Cu03 in 
Table [T| and the values computed earlier in Ref. |TT] We 
presume that these large differences, e.g., about 0.5 eV 
for the yz hole states, are mainly related to the approx- 
imations used in Ref. |TT] for representing the Cu and O 
ions on NN plaquettes. While here we represent those 
species at the all-electron level, in the earlier stud}/^^ the 
NN Cu and O ions were modeled by Mg^+ ions and for- 
mal 2— point charges, respectively. Point charges were 
also used for the embedding in Ref. [11] 

We note at this point that in an AF lattice the total 
energy of a given state within the d^ manifold is a sum 
of a crystal- field contribution, i.e., an on-site crystal- field 
splitting, and a magnetic term (see also the discussion in 
Refs. [HI and [TO]). As mentioned above, the AF NN spin 
coupling constant J is remarkably large in Sr2Cu03 and 
SrCu02, 0.20 to 0.25 eV^^ From the exact Bethe- 
ansatz solution for the ID Heisenberg Hamiltonian,l^SMl 
the AF ground-state stabilization energy is J In 2. On 
the other hand, from overlap considerations, we con- 
clude that for the crystal-field excited states the (su- 
per)exchange with the NN Cu dj.2_y2 spins is either zero 
or much weaker. For technical reasons, see the discussion 
in Sec. II, the quantum chemical calculations were here 
performed for a FM cluster. For a meaningful comparison 
between the MRCI results and experimental RIXS data, 
one should therefore subtract a term J In 2 from the rela- 
tive RIXS energies, representing the energy stabilization 
of the AF ground-state with respect to the crystal-field 
excited states. These considerations are relevant in light 
of future RIXS measurements on AF ID cuprates. 



B. Corner-sharing CuFe octahedra in KCuFa 

KCUF3 is a prototype material for systems with strong 
coupling among the charge, orbital, and spin degrees of 
freedom.'^^^ The crystalline structure of this compound 
is perovskite-like,'!^ i.e., three-dimensional (3D). The de- 
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TABLE II. Cu d-d excitation energies in KCuFs (eV). The 
ground-state Cu t2pC/^2c/^2_22 configuration is taken as refer- 
ence. The Jahn- Teller distortions occur within the xy plane. 
The MRCI values include Davidson corrections.^ 



Hole orbital 


CASSCF 


MRCI 


Experiment" 















0.76 


0.85 


0.71 - 1.02 


XZ 


0.89 


1.01 


1.05-1.15 


xy 


1.04 


1.17 


1.21 - 1.37 


yz 


1.11 


1.25 


1.34- 1.46 


Optical absorption, Ref. 43 


The numbers correspond to 



the onsets and the maxima of the absorption bands. 



generacy of the Cu Cg levels in an ideal perovskite lattice 
is lifted however through cooperative Jahn- Teller distor- 
tions. The latter imply a configuration with alternating, 
longer and shorter, Cu-F bonds for Cu ions along the x 
and y axes and orbital ordering in the xy plane. The hole 
in the Cu 3(i shell thus alternately occupies '^d^'^-z'^ and 
3(i^2_22 orbitals. 

In the xy plane, the magnetic couplings are weak and 
FM. On the other hand, along the z-axis the NN J is 
large and AF.^^ Actually this makes KCuFs a close to 
ideal ID Heisenberg antiferromagnet. The predictions 
for the lo w- lyin g spin excitations of the ID AF Heisen- 
berg chaii]^^^ have been confirmed by inelastic neutron 
scattering experiments."^ 

The magnetic behavior changes from ID to 3D at the 
Neel temperature Tn = 39 K. The emergence of sharp 
crystal-field absorption peaks in the optical spectra ap- 
proximately 10 K above Tn has been interpreted as ev- 
idence for a symmetry change related with a crossover 
from dynamic to static displacements of the F ions.*^ 
In this picture, the orbital and 3D AF ordering are in- 
timately related, in the sense that the former paves the 
road for the latter. 

The d-d excitation energies seen in the optical spec- 
tra were compared in Ref. 43 with the outcome of DFT 
band-structure calculations. However, strictly speaking, 
DFT is a ground-state theory. Here, we provide re- 
sults of excited-state CASSCF and MRCI calculations 
for the Cu level splittings in KCuFa, see Table V. The 
MRCI treatment brings corrections of 0.1-0.15 eV to the 
CASSCF splittings, somewhat smaller than for the ox- 
ides discussed above. This is related to the more ionic 
character of the Cu-F bond in KCuFa and smaller de- 
gree of Cu 3(i and F 2p orbital mixing. As a general 
trend, the MRCI results tend to slightly underestimate 
the values corresponding to the maxima of the experi- 
mental absorption peaks, by 0.1-0.2 eV. One obvious ef- 
fect of the less anisotropic environment is a much smaller 
splitting between the two eg components. This particu- 
lar electronic-structure parameter is in fact the smallest 
among the Cu d^ compounds investigated here. 



C. Chains of edge-sharing CUO4 plaquettes 

Cuprates in which the CUO4 plaquettes are edge- 
sharing, are characterized by weak superexchange inter- 
actions between NN Cu spins, which is due to Cu-O-Cu 
bond angles that are close to perpendicular. Depend- 
ing on the detailed crystal structure, the NN interac- 
tion is usu ally F M (LiVCu04, Li2Cu02)'^'^ or even AF 
(CuGeOs ) j ^'^ l ^^ l The in-chain next-NN exchange is always 
AF and causes frustration independently of the sign and 
size of the NN coupling. Further, the longer-range in- 
trachain and interchain couplings are sizeable as well in 
some of these compounds and lead to an intricate compe- 
tition of magnetic interactions. This has resulted in an 
appreciable scientific interest in the edge-sharing chain 
systems. 

In LiVCu04 the ratio of the FM NN coupling and the 
AF next-NN exchange - the frustration parameter - is 
controversial .^^^^ For LiCu202 even the sign of the NN 
coupling constant gave rise to debate.!^^^^ It has been 
also argued that a large fourth-nearest-neighbor AF cou- 
pling constant is responsible for the incommensurate heli- 
magnetic ground state of LiCu202!^In CuGeOs the NN 
exchange is AF and this system is famous to undergo a 
so-called spin-Peierls transition^^ to an AF state with 
a gapped excitation spectrum. When deriving the values 
of the shorter- and longer-range exchange interactions in 
terms of downfolding techniques, also the 3(i-level exci- 
tation energies come into play. 

The calculated CASSCF and MRCI results for the 
Cu d-level splittings in LiVCu04, CuGeOa, LiCu2 02, 
and Li2Cu02 are listed in Table [TTT) The finite embed- 
ded clusters on which the correlation treatment is car- 
ried out include one central and two NN CUO4 plaque- 
ttes for the single chain compounds LiVCu04, CuGeOa, 
and Li2Cu02. For the zigzag double-chain system 
LiCu2 02|^the cluster C includes four adjacent plaque- 
ttes. The coordinate framework is chosen such that the 




FIG. 2. Chain of edge-sharing Cu04 plaquettes in LiVCu04. 
The ground-state hole orbital is dxy in the chosen reference 
system. The six NN VO4 octahedra around the reference Cu 
site are highlighted in the figure. 
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TABLE III. CASSCF/MRCI d-d excitation energies for edge- 
sharing chains of Cu04 plaquettes (eV). The MRCI values 
include Davidson corrections.'^^ The ^/-axis is parallel to the 
Cu02 chains, z is perpendicular onto the Cu04 plaquettes, 
see text. 

Hole orbit al LiVCu04 CuGeOs LiCu202 Li2Cu02 

xy 

2 2 

X -y 
xz 



yz 

^2 





1.08/1.28 1.19/1.41 1.13/1.43 1.20/1.52 

1.25/1.47 1.42/1.61 1.58/1.88 1.72/2.04 

1.29/1.52 1.45/1.64 1.64/1.94 1.72/2.04 

0.98/1.18 1.50/1.71 1.67/1.98 1.93/2.30 



states employed in a fully ionic picture. When large devi- 
ations from the formal value of 2— occur for the effective 
charges of all adjacent O ions, the strongest affected is 
actually the ground-state energy because there are four 
O neighbors to which the lobes of the Cu dxy orbital are 
directed. This "destabilization" of the ground-state en- 
ergy of LiVCu04 explains the nearly uniform downward 
shift of the dx2_y2^ dyz^ and dxz excited hole states, by 
about 0.15 eV as compared to CuGeOs (see Table III). 



D. CuO 



y-ax.is is parallel to the Cu02 chains and z is perpen- 
dicular onto a given CUO4 plaquette. In other words, 
there is a rotation of 45° around z as compared to the 
standard coordinate framework for planar or octahedral 
coordination and the x and y axes are not along the in- 
plane Cu-0 bonds. The in-plane a-type Cu 3d orbital 
that is singly occupied in the ground-state configuration 
is thus denoted dxy. 

As for the corner-sharing plaquettes in Sr2Cu03 and 
SrCu02 (see Sec. III. A), inclusion of single and double 
excitations on top of the CASSCF wavefunctions brings 
a nearly uniform upward shift of 0.1 to 0.2 eV for the d- 
d transitions. By adding Davidson corrections, the rela- 
tive energies of the excited states further increase by 0.05 
to 0.1 eV. The MRCI values listed in Table ITITl include 
such Davidson correction terms. A second effect which 
deserves attention is the large variations found for the 
relative energy of the dz^ excited hole state. It is known 
that the length of the apical Cu-0 bond varies widely 
in Cu d^ oxides. The strong influence of the apical bond 
length on the d^^ hole excitation was previously stressed 
for 2D cuprates in, e.g., Refs. 9 and 10, Using simple 
electrostatic arguments, when the negative apical ions 
are closer to the Cu site, less energy is required for ex- 
citing the in-plane Cu 3d hole to the orbital pointing 
toward those apical ligands. As concerns the chain-like 
compounds from Table |III[ there are two apical O ions 
in CuGeOs, one apical in LiCu202, and no apical ligand 
in Li2Cu02. The relative energy of the d^^ hole state 
consequently increases from 1.71 eV in CuGeOs to 1.98 
in LiCu2 02 and to 2.30 eV in Li2Cu02. In LiVCu04 
there are two apical O ions as in CuGeOs but the apical 
bond lengths are much shorter as compared to CuGeOs, 
2.49 versus 2.76 A (see Refs. W and ITf). The d^y to 
dz2 excitation energy in LiVCu04 is thus substantially 
lower, 1.18 eV, even lower than in La2Cu04, where the 
apical bond length is 2.40 A and the transition to the 
dz2 orbital occurs at about 1.4 eV.^ However, one ad- 
ditional effect which comes into play in LiVCu04 is the 
strong covalency between the V ions and both in-plane 
and apical O neighbors of the Cu sites. Such covalency 
effects on the VO4 tetrahedra, see Fig. 2, give rise to 
large deviations from the formal V^^ and valence 



In spite of its apparent simple chemical formula cupric 
oxide, CuO, is a complex material. It was for instance 
only very recently discovered that it develops at the AF 
transition temperature Tn of 230 K also a net ferroelec- 
tric polarization, which is of substantial experimental^ 
and theoretical^^ interest. As the ferroelectricity is 
completely due to magnetism, CuO is classified as a type- 
II multiferroic.l^ The complexity of this material system 
is related to the fact that each of the distorted CuOe 
octahedra shares corners and edges with fourteen other, 
adjacent octahedra. As for most of the Cu oxides, on 
a given octahedron the two apical-like Cu-0 bonds are 
much longer than the other four Cu-0 bonds, 2.78 ver- 
sus 1.95 or 1.96 A.^^The quantum chemical calculations 
were therefore carried out on a cluster whose active re- 
gion Ca is defined by one reference CUO4 plaquette and 
the nearest fourteen Cu neighbors, see Fig. 3. The buffer 
region Cb includes the apical O ions of the reference Cu 
site and for each adjacent Cu ion the nearest four ligands. 




FIG. 3. Crystal structure of CuO and sketch of the cluster 
used for the calculations. The reference Cu and the fourteen 
adjacent Cu sites are shown in dark and light blue, respec- 
tively. O ions of the central octahedron are depicted in or- 
ange. The thin dashed line parallel to the y-axis connects 
the central Cu site and two Cu NN's along the same chain of 
edge-sharing CUO4 plaquettes. For each apical O, there are 
two NN Cu ions on adjacent chains parallel to the y-axis. 
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TABLE IV. CASSCF and MRCI d-d excitation energies 
for CuO, see text. The MRCI values include Davidson 
corrections.^^ All numbers are in eV. 



Hole orbital 


CASSCF 


MRCI-1 


MRCI-2 


xy 











2 2 

X -y 


1.13 


1.37 


1.38 


yz 


1.54 


1.75 


1.79 


xz 


1.60 


1.81 


1.85 


z' 


1.71 


1.93 


1.96 



CASSCF and MRCI results for the Cu d-level splittings 



in CuO are given in Table [IV| CASSCF calculations were 
first performed for a FM configuration of the fifteen Cu 
sites. A non-standard local reference system was chosen 
for the central Cu site, with a rotation of 45° around z as 
for the cuprates discussed in Sec. III.C. The ground-state 
hole orbital is dxy in this framework. CASSCF relative 
energies for the crystal-field excited states are listed in 
the second column of Table HVl 

In a next step, MRCI calculations were carried out. 
Due to the large number of open-shell NN Cu sites, for 
CuO we were forced to restrict the CAS orbital reference 
space to the set of 3d orbitals of the central Cu ion. Two 
different types of approximations were employed. In a 
first set of restricted multiconfiguration calculations, out 
of the fourteen linear combinations of NN Cu dxy func- 
tions, we froze the occupation of the seven lowest-energy 
(bonding-like) orbitals to 2 while the occupation of the 
higher-lying (antibonding-like) orbitals was restricted to 
0. Only the Cu 35, 3p, 3d and O 2p electrons on the cen- 
tral Cu04 plaquette were correlated in the subsequent 
MRCI treatment. The difference between the CASSCF 
and MRCI values for each particular state in this set of 
CASSCF and MRCI calculations was added afterwards 
to the initial CASSCF d-d splitting. The resulting num- 
bers are listed in the third column of Table IIVI In the 
second set of MRCI calculations, we replaced the Cu^+ 
d^ neighbors by closed-shell Zn^+ d^^ ions. Those MRCI 
correlation induced corrections were again added to the 
initial CASSCF d-d excitation energies and the corre- 
sponding results are given in the fourth column of Table 
IV As shown in Table |IV[ the MRCI treatment brings 
shifts of 0.2 to 0.25 eV to the CASSCF splittings. In the 
two approximation schemes, the MRCI corrections for 
each particular state are nearly the same, i.e., the MRCI 
correlation energies do not depend on details of the NN 
3d electron configuration. 

The lowest crystal- field excitation is from dxy to dx2-y2 
and estimated in Table [TVl at about 1.4 eV. This value 
compares well with the measured optical gap, 1.3-1.6 
Q-yMMMi RIXS measurements on CuO by Ghiringhelli et 
a/P^ show that most of the weight of the d-d excitation 
spectrum is between 1.7 and 2.3 eV, with clear peaks at 
1.85 and 2.15 eV. Our MRCI data seem to somewhat un- 
derestimate those RIXS results, even when a J In 2 cor- 



rection is considered to account for the superexchange 
interactions along the AF chains in CuO.^^ More flexible 
Cu and O basis sets are expected to further reduce these 
small differences between theory and experiment. 



E. A ladder cuprate: CaCu2 03 

In the so-called ladder compounds, the transition- 
metal ions are arranged in a planar network of ladders 
which can display two or more legs. These systems 
became subject of intense study when Dagotto et al^ 
found theoretical evidence that the isolated 6' = l/2 AF 
two- leg ladder has a finite spin gap, i.e., a finite energy 
is needed to create a spin excitation. 

Physical realizations are found in the vanadium oxide 
compound CaV205^^ and in cuprates such as CaCu203, 
SrCusOs,^ and (Sr,Ca)i4Cu2404i.^ The latter com- 
pound also displays superconductivity.'^ 

We analyze in this section the d-level electronic struc- 
ture of CaCu203. For most of the two- leg ladder 
cuprates, including CaCu203, octahedra on the same lad- 
der share common corners while adjacent octahedra on 
different ladders share common edges. One peculiar fea- 
ture of CaCu203 is that on each rung of a ladder the 
Cu-O-Cu angle is 123°, much smaller than in other lad- 
der compounds.'^ This is the reason the intersite mag- 
netic couplings across the rungs are AF and small, about 
10 meV, while the superexchange along the leg of the 
ladder is AF and one order of magnitude larger, about 
135 meVW^ In other words, from the magnetic point of 
view the system is rather ID, with weakly interacting 
AF Cu02 chains. Choosing the a-axis parallel to the 
rung and b along the leg of the ladder, the ground-state 
hole orbital is dx'2-y2. 

The quantum chemical calculations were carried out 
on a cluster whose active region Ca is defined by one 
reference CUO4 plaquette and the five adjacent Cu ions 
(one Cu NN on the same rung, two leg NN's, and two 
Cu NN's on an adjacent ladder). The buffer region Cb 
includes the two apical O ions of the reference Cu site, 
the nearest four ligands for each adjacent Cu ion, and six 
Ca neighbors. The apical O ligands are at rather large 
distance from the central Cu site, 3.0 A.^ 

CASSCF and MRCI results for the Cu (i-level splittings 



TABLE V. CASSCF and MRCI d-d excitation energies for 
the ladder system CaCu2 03 (eV). The MRCI values include 
Davidson corrections.^^ 
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FIG. 4. MRCI Cu d-d excitation energies in ID Cu d^ compounds. Depending on the reference coordinate system, see text, 
the symmetry of the lowest t2g hole state (in red) is denoted either as dxy (corner-sharing chains of plaquettes or octahedra) 
or d^2_y2 (edge-sharing chains). For KCuFs, the Cg excited state (in orange) has symmetry. For all other compounds, the 
Cg excited state has symmetry. 



in CaCu203 are listed in Table [Vl The MRCI treatment 
and the Davidson corrections bring shifts of 0.2 to 0.3 
eV to the CASSCF d-d excitation energies. The excita- 
tion energies to the xz, yz, and z^ orbitals are somewhat 
smaller than in Sr2Cu03, SrCu02, and Li2Cu02 because 
in the latter compounds there are no apical ligands. At 
the same time, the splitting between dx2_y2 and is 
larger than in the other cuprates investigated here, see 
Fig. 4, because the apical Cu-0 bond length is the largest 
in CaCu203. 



IV. CONCLUSIONS 

It has been shown here that wavefunction-based 
electronic-structure calculations are well-suited for com- 
puting local charge excitations as observed in transition- 
metal oxides by optical and RIXS measurements. While 
this was demonstrated previously for a number of lay- 
ered cuprates,^ we have extended here those applica- 
tions to ID and ladder compounds. Results and trends 
for the d-d excitations in low-dimensional cuprates are 
summarized in Fig. 4. The ab initio data compare well 
with the results of optical absorption measurements on 
KCuFsF^and somewhat underestimate the excitation en- 
ergies found by RIXS experiments for CuO.^^ Future 
RIXS experiments on other ID cuprates included in 
Fig. 4 will constitute a direct test for the symmetry and 
energy of the different d-d excitations that we predict 
here. 



From a more general perspective, the results reported 
here clearly indicate the significant potential of the ab 
initio wavefunction-based methods in the context of cor- 
related electronic materials and the need to further de- 
velop such techniques in parallel with approaches based 
on density functional theory. The latter has revolution- 
ized the field of electronic-structure calculations but its 
limits are also clear. Computations of excited states of 
strongly correlated electrons are in a way the worst case 
for DFT. One is then forced to avail oneself to alter- 
natives, i.e., wavefunction-based methods. An attrac- 
tive feature is that strong correlations can be handled 
by CASSCF calculations even when we deal with infinite 
periodic systems. In all approximations which are being 
made, the accuracy of the wavefunction-based methods 
can be progressively increased. Here we supplemented 
the CASSCF calculations by MRCI, yet, alternative sup- 
plementary methods are possible. Further exploration of 
such ab initio techniques is planned for the future. 
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